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Abstract 

General relativity is a covariant theory of two transverse, traceless graviton degrees 
of freedom. According to a theorem of Hojman, Kuchaf , and Teitelboim, modifications 
of general relativity must either introduce new degrees of freedom or violate the prin- 
ciple of general covariance. In this paper, we explore modifications of general relativity 
that retain the same number of gravitational degrees of freedom, and therefore explic- 
itly break general covariance. Motivated by cosmology, the modifications of interest 
maintain spatial covariance. Demanding consistency of the theory forces the physical 
Hamiltonian density to obey an analogue of the renormalization group equation. In 
this context, the equation encodes the invariance of the theory under flow through 
the space of conformally equivalent spatial metrics. This paper is dedicated to set- 
ting up the formalism of our approach and applying it to a realistic class of theories. 
Forthcoming work will apply the formalism more generally. 
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1 Introduction 



For nearly a century, general relativity has been the most successful paradigm for interpreting 
and understanding classical gravitational phenomena, and to this day there have been no 
unequivocal refutations of general relativity Nonetheless, there are compelling reasons to 
study alternative gravitational theories. 

Perhaps the most obvious reason is to explain empirical anomalies, most notably the 
observed cosmic acceleration. While it is true that this phenomenon can be understood 
in terms of a cosmological constant in pure general relativity, this approach has the draw- 
back that we have no handle on what physics might set the magnitude of the cosmological 
constant. It is currently an outstanding theoretical challenge to determine what physical 
degrees of freedom are associated with late-time acceleration. Dynamical theories of dark 
energy postulate scalar degrees of freedom similar to those invoked to account for inflation. 
Unfortunately, there is no unambiguous evidence for such scalar degrees of freedom. One mo- 
tivation for this paper is the possibility that cosmic acceleration might be directly associated 
with the graviton degrees of freedom. 

Apart from any attempt to understand empirical anomalies, there remains a compelling 
theoretical reason to study alternatives: to determine which features of general relativity are 
essential to its experimental success, and which features are merely incidental. To analyze 
the theory in this manner, we must know what freedom we have to modify the theory while 
retaining its explanatory power. 

The two transverse, traceless graviton degrees of freedom are a key feature of general 
relativity. Though graviton exchange has never been measured and gravitational waves have 
never been directly detected, there is substantial indirect evidence for the existence of these 
degrees of freedom. It is logically possible that additional gravitational degrees of freedom 
exist, but at the same time there is no unambiguous evidence for them. It is therefore natural 
to ask whether and how we can modify general relativity while preserving the same number 
of degrees of freedom. 

In this paper, we construct manifestly consistent modified theories of gravity that retain 
the same local degrees of freedom as general relativity. To evade the consequences of the 
theorem that general relativity is the unique theory of a massless spin-2 particle [1, 2], our 
theories break local Lorentz symmetry explicitly. Theories in which Lorentz symmetry is 
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only broken spontaneously necessarily rely on additional local degrees of freedom. These 
appear in the broken phase as massless Goldstone modes; an example of such a theory is 
ghost condensation [3]. 

Perhaps the biggest obstacle to modifying general relativity is the opaqueness of the 
theory. General relativity as formulated by Einstein and Hilbert is a covariant theory, which 
means that the equations of motion for the spacetime metric g^ u take the same form in any 
coordinate system. Unfortunately, invariance under coordinate transformations implies that 
the theory contains a great deal of gauge arbitrariness, and the underlying dynamical degrees 
of freedom of the theory have proven difficult to isolate. In fact, the notorious elusiveness of 
the physical degrees of freedom is partially responsible for the difficulty of quantizing general 
relativity. 

This gauge arbitrariness can be understood most clearly by treating general relativity 
as a constrained field theory. By writing the metric g^ u in ADM form 1 and discarding a 
boundary term, the Einstein-Hilbert action can be rewritten in canonical form as a theory of a 
spatial metric and a conjugate momentum tensor n i: * subject to four first class constraints 
T-Lfj, [4]. While this form of the theory is not written in terms of manifestly diffeomorphism 
covariant objects, the covariance of the theory can still be inferred from the covariance 
"algebra" satisfied by the which act as the generators of spacetime diffeomorphisms 
under the action of the Poisson bracket [5] . By representing the gauge symmetries of general 
relativity as constraints on the phase space of the theory, it becomes straightforward to 
count degrees of freedom. According to the standard counting prescription, it follows from 
the presence of four first class constraints in a theory of six canonical coordinates that 
general relativity contains two local degrees of freedom. 2 In the passage to quantum theory, 
these transverse, traceless degrees of freedom become the two polarizations of the graviton. 

To isolate the physical graviton degrees of freedom, one would have to solve the four 
constraints By taking the configuration space for the metric to be Wheeler's superspace, 
it is possible to solve the three momentum constraints %i by fiat, but the Hamiltonian 
constraint Hq has thus far defied solution in general. Unless the Hamiltonian constraint can 
be solved, the gauge-arbitrariness of general relativity cannot be eliminated. Fortunately, 
though no general solution to the Hamiltonian constraint has been found, it is possible to 



1 i.e., in terms of a spatial metric hij, a lapse N = N°, and a shift N % . 
2 See section 2.1 for more detail. 
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solve it in certain circumstances by imposing an appropriate gauge-fixing condition. 

In 1974, Hojman, Kuchaf, and Teitelboim (HKT) proved that general relativity with a 
cosmo logical constant is the unique covariant theory of a spatial metric and its conjugate 
momentum 7r* J ; in fact, general relativity is the minimal representation of the covariance 
algebra [6, 7]. It follows that alternative covariant theories of gravitation must introduce 
additional degrees of freedom beyond the two graviton degrees of freedom in general relativ- 
ity [8]. Conversely, alternative theories of a spatial metric and its conjugate momentum 
7r u cannot be covariant. To modify general relativity, one must either introduce new degrees 
of freedom or violate the principle of general covariance. 

We wish our theories to retain the same local degrees of freedom as general relativity, 
so in accordance with the theorem of HKT, our theories cannot be diffeomorphism covariant. 
This aspect of our approach is not necessarily a defect. The covariance property of general 
relativity may be intellectually appealing, but the form invariance of the equations of motion 
is purchased at the cost of substantial gauge redundancy. Moreover, since we do not observe 
exact spacetime symmetry in our universe, this property of general relativity is not neces- 
sarily key to the success of the theory. Simply put, on cosmological scales there is a strong 
asymmetry between the past and the future, and the observable universe has a preferred rest 
frame; these observations are conventionally understood as a result of spontaneous symmetry 
breaking, but explicit symmetry breaking is another logical possibility. 

That being said, on cosmological scales in the cosmological rest frame there is substan- 
tial evidence for spatial homogeneity and isotropy. To maximize the verisimilitude of our 
treatment, the theories we consider will retain explicit covariance under spatial diffeomor- 
phisms. To summarize, we will attempt to modify general relativity while preserving 1) the 
number of graviton degrees of freedom, and 2) covariance under spatial diffeomorphisms. In 
this paper, we develop a general framework within which to explore the freedom we have to 
modify general relativity while retaining these two desirable properties. 

Concretely, we will begin by recasting general relativity in spatially covariant form, 
by solving the Hamiltonian constraint while preserving the momentum constraints. We will 
solve the Hamiltonian constraint by imposing a cosmologically motivated gauge constraint: 
we will take the determinant of the spatial metric to be the measure of time. This operation 
destroys the manifest diffeomorphism covariance and local Lorentz covariance of the theory. 
We emphasize that this gauge breaks down in the general case when the determinant of 
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the spatial metric is allowed to evolve non-monotonically, but it is a natural choice when 
considering perturbative corrections to FRW spacetime. By solving the Hamiltonian con- 
straint, the determinant of the spatial metric and the trace of the momentum tensor drop 
out of the phase space of the theory. We thereby obtain general relativity as a theory of a 
unit-determinant metric and a traceless conjugate momentum tensor n 1 ^ subject to three 
first class momentum constraints "Hj, which act as the generators of spatial diffeomorphisms. 
By the standard counting prescription, the presence of three first class constraints T-Li in a 
theory of five canonical coordinates guarantees that spatially covariant general relativity 
contains two degrees of freedom, as it should. 3 

Our strategy for modifying general relativity relies on the fact that any theory of 
five canonical coordinates subject to three first class constraints contains two degrees of 
freedom. To modify general relativity, we will modify the functional form of the physical 
Hamiltonian density on the reduced phase space (hij, 7r u ), subject to the condition that the 
momentum constraints "Hj remain first class; to ensure the consistency of the modification, 
we will also demand that the constraints "Hj remain preserved by the equations of motion. 
Any theory that satisfies these two restrictions will retain manifest spatial covariance, and 
by the counting prescription will necessarily contain two graviton degrees of freedom. In 
this paper, we introduce the formalism necessary to pursue this program of modification 
and apply the formalism to a class of realistic theories. Forthcoming work will apply the 
formalism developed here to a broader class of theories [9]. 

The literature abounds with many and varied approaches to the pursuit of modified 
gravity theories, but covariant modifications of general relativity that introduce additional 
degrees of freedom have been the most widely explored. The well-known method for find- 
ing covariant theories is to construct a diffeomorphism-invariant Lagrangian density out of 
manifestly covariant objects by contracting all free spacetime indices. Using this technique, 
all manner of theories have been explored: scalar-tensor theories [10], theories with higher- 
order curvature terms [11, 12, 13], theories of massive gravity [14, 15, 16, 17, 18, 19], higher- 
dimensional gravity theories [20, 21, 22, 23], galileons [24, 25, 26], chameleons [27, 28, 29, 30], 
symmetrons [31, 32], etc. For a comprehensive review of Lorentz-invariant massive gravity 
theories with detailed references, see [2]. For a comprehensive review of observational tests 
of modified gravity, see [33]. 

3 See section 3.4 for more detail. 
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Non-covariant approaches have been tried as well, but there is no single unifying proce- 
dure for the construction of such theories. In general, the natural procedure for understand- 
ing non-covariant theories depends on which symmetries survive in the theory. For example, 
in [34] Lorentz-violating massive graviton theories were classified by assuming the graviton 
mass to be invariant under the three-dimensional Euclidean group. A prominent recent 
example of a non-covariant metric theory is Hofava-Lifshitz gravity [35, 36, 37] , 4 in which 
the phase space constraints are chosen to satisfy a non-relativistic version of the covariance 
algebra. Also of note is the work of Barbour and collaborators on theories of conformally 
equivalent spatial metrics [38]. 

This paper is organized as follows. In section 2, we cover the basic concepts of con- 
strained field theory in the context of analyzing the phase space and constraint structure of 
general relativity. In section 3, we show how to impose our cosmological gauge condition 
and solve the Hamiltonian constraint to obtain a consistent spatially covariant formulation 
of general relativity. In section 4, we introduce the formalism of our approach to modifying 
gravity in the context of theories with an ultralocal physical Hamiltonian density. In sec- 
tion 5, we apply our method to derive consistency relations for a class of realistic theories 
which includes general relativity. 



2 General Relativity as a Constrained Field Theory 



In this section, we will analyze general relativity by treating it as a constrained field theory. 
In particular, we will examine its phase space and constraint structure, and count its local 
degrees of freedom. 

Our starting point is the Einstein-Hilbert action with a cosmological constant, 



From this action, the general covariance of the theory is manifest, but the counting of degrees 
of freedom is not. The metric tensor has ten components, but the theory has only two 
independent local degrees of freedom. To facilitate the counting of degrees of freedom, it is 
conceptually simplest to rewrite the action in a manner which makes the counting manifest, 

4 Thc original incarnation [35] of Hof ava-Lifshitz gravity struggled with consistency issues [36] which were 
resolved in [37] by imposing a consistent constraint algebra. 




(1) 
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i.e., canonical form. To this end, the metric g^ v must first be expressed in ADM form, in 
terms of a lapse N = iV°, a shift N\ and a spatial metric Kif 

ds 2 = g^dx^dx" = -N 2 dt 2 + h l3 (dx l + JVMt) (dx j + iV%) . (2) 

Up to a boundary term, the action of general relativity is 

S = J dt d 3 x VhN (K tJ K ij - K 2 + R - 2 A) . (3) 

In this expression, indices are lowered with hij and raised with its inverse h ij , R = R {3) is 
the Ricci scalar of the metric h^, the extrinsic curvature tensor is defined by 

K* = \ N ' 1 ( h iJ - ViNj - VjN t } , (4) 

K = h lj Kij, and V, = V- 3 ' 1 is the covariant spatial derivative with respect to the metric h^. 
To obtain the canonical action, one must first define the momentum conjugate to the spatial 
metric 

n ij = ^ = Vh (K ij - Kh ij ) ■ (5) 
Shij 

the momentum is three-tensor density of unit weight. 5 By inverting the relation between 
and K % \ one can rewrite the action of general relativity in canonical form as 

S = J dtd 3 x (n^hij - N^U^j , (6) 

where 

Ho = -Vh(R - 2A) + i= (n^ - \^\f^j , 
U t = -2h tJ V k ^ k . (7) 
Variation of the action with respect to and yields Hamilton's equations, 



where the Hamiltonian H is 

H 



= J d 3 xN^. (9) 



5 According to the standard convention, the weight of a tensor density is the number of times \fh multiplies 
the underlying tensor. 
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To evaluate the above variational derivatives, one must use the relations 



Shjjjx) 
Shki(y) 



5rr kl (x) 



6%6 3 (x-y), 



(10) 



where 




2 m + • 



(ii) 



Defining the Poisson bracket 




Sh mn (z) Sn mn (z) Sn mn (z) Sh mn (z) 



5A 6B SA SB 



(12) 



the equation of motion for any quantity A(hij, ir^ , t) can be written as 



A 



dA 

~dt 
OA 

~dt 



+ {A,H} 




(13) 



If A has no explicit dependence on time, its evolution is generated by its Poisson bracket 
with the Hp. 

Variation of the action with respect to iV yields the four constraints 



The symbol ~ denotes weak equality, or equality after the constraints % M ~ have been 
enforced. For example, if X = Y + X^V,^, then X ~ Y. Since the constraints define a surface 
in phase space, weak equality is also termed equality on the constraint surface. As an aside, 
it follows from (9) and (14) that H ~ 0; the vanishing of the Hamiltonian on the constraint 
surface is a feature common to covariant theories. 

There is no tt^N 11 term that would allow us to compute a variational expression for 
iV M , so the time evolution of iV M is unconstrained by the action. The four functions iV M are 
thus arbitrary until and unless we gauge-fix them. 

2.1 Constraint Properties & Degrees of Freedom 

Before examining the constraints more closely, we pause to review some terminology first 
introduced by Dirac for describing constrained theories. A quantity whose Poisson bracket 



~ 0. 



(14) 



7 



with each of the constraints vanishes (identically or weakly) is termed first class; a quantity 
whose Poisson bracket fails to vanish weakly with at least one constraint is termed second 
class. A first class constraint has vanishing Poisson bracket with all constraints, while a 
second class constraint has non- vanishing Poisson bracket with at least one other constraint. 
In most cases of interest, first class constraints generate gauge symmetries under the action 
of the Poisson bracket. Second class constraints can usually be solved, either implicitly (by 
using the "Dirac bracket") or explicitly (by expressing some phase space variables in terms 
of others). 

In general relativity, the constraints generate spacetime diffeomorphisms. By direct 
calculation — see appendix A for details — it is possible to prove that the constraints % M 
are first class, {H^x), 1-L u (i))} ~ 0. This means that the symmetry generators are closed 
under the action of the Poisson bracket, as they should be in order to consistently represent 
any kind of a symmetry. In particular, 

{Ho(x), U (y)} = W(x)d xi 5 3 (x - y) - W(y)d y ,5 3 (x - y) , 
{Ho{x),Hi{y)} = U (y)d x .6 3 (x - y) , 

{Hi(x),Hj(y)} = nj(x)d xi 5 3 (x - y) - H t (y)d y] 5 3 (x - y) . (15) 

Obeying this covariance "algebra" is the necessary and sufficient condition for the to 
consistently generate spacetime diffeomorphisms [5]. In fact, four first class constraints 
obeying this algebra are guaranteed to arise in any covariant field theory; in canonical form, 
general covariance of the action is encoded in this constraint algebra. 

For the constraints to be consistent with the equations of motion, the constraints must 
be preserved by the equations of motion, i.e., "H M ~ 0. Since &H^/dt = 0, applying the 
equations of motion to T-L^ yields 

n,(x) = J d 3 yN»( y ){H»(x),H v (y)}- (16) 

From the first class character of the constraints, it follows that Iri^ ~ 0, as desired. 

The Hamiltonian formulation of GR is a theory of a spatial metric hij and its conjugate 
momentum so the theory contains twelve canonical (or six real) variables. However, 
these variables are not independent. First, they are related by the four constraints ~ 0. 
Second, the equations of motion for and ir^ depend on the four arbitrary functions N^ 1 ; 
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to gauge-fix the iV M would require imposing four gauge-fixing constraints. 

6 • hij 's + Q-ir ij 's -A-H/s - 4 • N»'s = 4 canonical DoF . (17) 
The theory therefore has four canonical (or two real) degrees of freedom. 



3 Spatially Covariant General Relativity 

We would like to depart from general relativity by modifying the equations of motion for the 
two graviton degrees of freedom. Ideally, we would like to solve all four gauge constraints, 
go down to the physical phase space, and modify the theory at that level. In this way, 
we would circumvent all the difficulties of consistently modifying a constrained field theory. 
Unfortunately, we do not know how to do this. 

One possible approach is to modify the equations of motion for the phase space variables 
hij and n 1 *. However, the counting of degrees of freedom in general relativity relies on the 
fact that the four constraints "H M satisfy a consistent first class algebra, namely the covariance 
algebra of equation (15), and we know from the HKT theorem that any modification of the 
action for and n lj will destroy this algebra. If we modify the action for the phase space 
variables and , we must impose an alternative constraint structure that consistently 
constrains the phase space to the same degree as the covariance algebra; this is the approach 
taken in [37]. Since we take the point of view that full spacetime covariance is a spurious 
symmetry, we do not wish our theory to contain a constraint structure that implies the same 
degree of redundancy as the covariance algebra. 

Though spacetime symmetry is manifestly broken on cosmological scales (whether 
spontaneously or explicitly), there is strong evidence for spatial homogeneity and isotropy, so 
we will attempt to modify general relativity while preserving the manifest spatial covariance 
of the theory. To obtain a spatially covariant formulation of general relativity to modify, we 
will solve the Hamiltonian constraint "Ho while leaving the three momentum constraints Hi 
intact. The Hamiltonian constraint is famously hard to solve in general, but we are inter- 
ested in using our theories in a cosmological context, so we will solve it using a gauge-fixing 
constraint which is well-defined on an expanding FRW background. 
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3.1 Metric Decomposition 

Before gauge-fixing, we decompose the metric hij into a conformal factor = h 1 ! 3 and a 
unit-determinant metric h^, i.e., 

= Qhij . (18) 

Note that f2 = (\fh) 2 / 3 is a three-scalar density of weight 2/3, while h^ is a three-tensor 
density of weight —2/3. The scalar density we will work with is not the conformal factor Vt, 
but the volume factor u = \fh = f2 3 / 2 , which is a scalar density of unit weight. We choose 
ui because its conjugate momentum, 

is a three-scalar and hence invariant under a rescaling of Q or u; this fact will simplify 
matters in sections 4 and 5. The momentum conjugate to is 

~ij = j^L = Q ^ij _ hi l3 it k ^j = tutt (^K ij - ^Kh 1 ^ , (20) 



which is a traceless three-tensor density of weight 5/3; the quantity 



7T 



is the corresponding traceless three-tensor. By defining the traceless projection tensor 5™ 

= ^-\h l] h k \ (22) 



we can write n 1 ^ more compactly as 

r j = tt5%n kl = ttu~5%K kl . (23) 
The phase space variables and n lj can thus be written as 

hij = u 2/ %, ir ij = cu~ 2/3 r j + ^tiiu 1/3 TT u . (24) 

The decomposition of the spatial metric into a volume factor and a unit-determinant met- 
ric is completely general. Though the corresponding conjugate momenta were derived by 
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taking variational derivatives of the Einstein-Hilbert Lagrangian, the decomposition of the 
momentum tensor into its trace part and its traceless part is likewise completely general. 
Those familiar with the techniques of numerical relativity may be reminded of the York- 
Lichnerowicz conformal decomposition or the BSSNOK (Baumgarte, Shapiro, Shibata, Naka- 
mura, Oohara, and Kojima) formalism [39]. 

3.2 Cosmological gauge 

To solve the constraint Ho, we must first gauge- fix the lapse N with a gauge- fixing constraint 
X for which {Ho, x} 00 0; this renders Ho second class, and hence solvable. Since we wish to 
retain explicit spatial covariance, our constraints Hi must remain first class. 

In a cosmological context, it is natural to use the volume factor of the spatial metric 
as a clock, so that t = t(u); we call this cosmological gauge. As mentioned in the intro- 
duction, cosmological gauge is only valid when the determinant of the spatial metric evolves 
monotonically, so this procedure is only valid when considering perturbative corrections to 
FRW spacetime. When the evolution of u is monotonic, t(u) is an invertible function, so 
this gauge is equivalent to taking the volume factor u; to be a function of time, i.e., u> = u(t). 

To impose cosmological gauge, we add to the canonical action of general relativity a 
gauge-fixing constraint 

X = uj - u(t) , (25) 

along with a corresponding Lagrange multiplier A. The new gauge-fixed action is 

S' = J dtd 3 x (n ij h i:j - N^H? - Ax) . (26) 

Varying the action with respect to A then reproduces the constraint 

X~0. (27) 

By direct calculation — see appendix B for details — one can verify that 

{Ho(x), X (y)} = ln\(x)5 3 (x~y); (28) 

the constraints Ho and x are thus second class, so we expect to be able to solve them. The 
only wrinkle is that 

{Hi(x), X (y)} = Vh&)d xi 6 3 (x - y) , (29) 
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so the constraints Hi are also second class! By shuffling our constraints slightly, we can 
obtain a set of two second-class constraints and three first class constraints, and thereby 
preserve explicit spatial covariance. Indeed, since 

fey/htxjd* (^|y) Mv)} ~ VW)d xi 5 3 (x-y), (30) 
it follows that the combination 

n i = n l -2Vhd l 

= % i -2VhV i (^f ) j (31) 

obeys 

[n^Sjiv)} ~ {Hi(x),Ho(y)} ~ {^(x), X (y)} ~ 0. (32) 

The interpretation of this result is simple. The HiS generate spatial diffeomorphisms, while 
Ho generates time translation. A generic spatial diffeomorphism will alter the conformal 
factor of the spatial metric. If the conformal factor is taken to be the measure of time, 
then the Hi's, by altering the conformal factor, will generate time translation, while Ho, 
by generating time translation, will alter the conformal factor. The "Hj's generate spatial 
diffeomorphisms that preserve the conformal factor, so they must differ from the H^s by the 
gradient of a compensating time translation term. 

From the definition of Hi, it is apparent that demanding x ~ and H^ ~ is equivalent 
to demanding x ~ 0, Ho ~ 0, and Hi ~ 0. The latter set of constraints has the virtue that 
the Hi are first class, and can thus consistently represent symmetries. We therefore take 
our five constraints to be the two second class constraints \ an d Ho and the three first class 
constraints Hi- Using Hi — Hi + 2^/hVi (Ho/^k), the gauge-fixed action can be rewritten 
in terms of Hi as 

S' = J dtd 3 x (n^hij - N°H - N l Hi - 2VhN i V i - X^j . (33) 

Upon integration by parts, the action becomes 

S'= I dt d 3 x (ix i3 h l3 - NHq - N { Hi - Ax) , (34) 
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where 



\fh 

N = N — 2^rViN l . (35) 



71 k 



Variation of the action S' with respect to and 7r y yields Hamilton's equations, 

5H' 5H' 

k - {x) = 5^(x)> r3{x) = -Jh-(x)> (36) 



where the new Hamiltonian H' is 



H' = J d 3 x (nH + N'Ui + Ax) . (37) 



The equation of motion for any quantity A(hij,ir l i,t) is therefore 



OA 
~dt 



+ J d 3 y (N{y){A,U {y)} + N\y){AMy)} + Ky){^x{y)}), (38) 



where the Poisson bracket is defined as in (12). Variation of the action S' with respect to 
N, A, and N l yields the five constraints 

^ ~0, x~0, Hi~0. (39) 

The action does not contain time derivatives of the Lagrange multipliers, so at first their 
evolution appears unconstrained. Since the %i are first class, the three functions N l are 
indeed arbitrary until and unless we gauge-fix them. The evolution of N and A, however, 
will be determined by demanding the consistency of H and x with the equations of motion. 

For the constraints to be consistent with the equations of motion, they must be pre- 
served by the equations of motion; we therefore demand that %i ~ 0, "Ho ~ 0, and x ~ 0. 
Since the Irii are first class and d{Ki = 0, it follows at once that "Hj ~ 0. Since dtHo = 0> 
{Ho(x), Hoiy)} ~ 0, and {Ho(x),'Hi(y)} ~ 0, it follows that 

Ho(x)~ J d 3 y\(y){H (x),x(y)} 

~ \x(x)7r\(x) . (40) 

On a flat FRW background, 6 K = 3d/ a and hence 7r\ = —2uK = —6aa 2 . Since we are only 
considering gravity on an expanding background, we assume that tt\(x) oo more generally. 



6 A spatially-flat FRW spacetime corresponds to N = 1, Ni = 0, and hij — a 2 (t)Sij. 
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The demand Ho ~ thus implies 



A ~ . 



(41) 



Since {x{x),x{y)} = 0, {x(x) ,Ui(y)} ~ 0, and dx/dt 



■oj(t), it follows that 



X(x)~-u>(t) + J d 3 yN(y){x(x),n (y)}, 
~-dj(t)-\N(x)n\(x). 



(42) 



Since n l i oo 0, demanding x 



allows us to solve for N, 



N ~ 



2w(t) 



(43) 



The functions N and A are thus not arbitrary. Since N = N + 2Vh(V i N i )/n k k , the lapse N 
has not been completely gauge-fixed, but its arbitrariness stems solely from its dependence 
on the three arbitrary functions N l . 

As a check, let us revisit the counting of degrees of freedom in cosmological gauge. For 
these purposes, the only effect of gauge-fixing is to replace the first class constraint Ho ~ 
and the arbitrary function N with the second class constraints Ho ~ and x ~ 0. This 
modifies the left hand side of equation (17), but does not change the final tally. 

6 • h i:j 's + 6 • 7r ij 's - 1 • Ho - 1 • x - 3 • Hi 's - 3 • N* 's = 4 canonical DoF . (44) 

After gauge-fixing, the theory still has four canonical (or two real) degrees of freedom. 

3.3 Solving Ho and x 

In this section, we will solve the constraints T-Lq and x t° obtain a spatially covariant for- 
mulation of general relativity as a theory of a unit-determinant metric and its conjugate 
momentum ti % K This will set the stage for modifying general relativity in section 4. 

Since x an d "Ho ar e second class, they can be solved explicitly to yield expressions for 
uj and 7r w in terms of t, hij, 7f u , and spatial derivatives. "Solving" for u is trivial: u = u(t). 
Solving for ix w requires us to take a square root and pick a sign, which amounts to picking 
either an expanding or a contracting background. We pause to emphasize once again that 
our procedure is only valid in a cosmological context, when the conformal factor of the spatial 
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metric can be assumed to be evolving monotonically. To pick the sign corresponding to an 
expanding background, first recall that 

n u = ~K. (45) 

On a flat FRW background, K = 3a/ a and hence tt^ = —4a/ a. An expanding FRW 
background therefore corresponds to 7r w < 0. Returning to the general case, we choose 
7r w < to obtain 



I „ = TGRS -,/-y^-- 375 + 2A, (46) 

where indices are raised and lowered with h^, and R is the Ricci scalar for h^. Substituting 
these results for u and back into the action S' yields the action of general relativity on 
the reduced phase space (hij, 7r u ), 

S" = J dt d 3 x ( Tt ij hij + n^cu - ITili ) , (47) 

where 

Hi = -2hijV k n jk - uVi^ , (48) 
and Vj is the covariant derivative with respect to hij. This action yields the new Hamiltonian 

H" = J d 3 x ( -cun^ + N l H t ) . (49) 

The term n^hij has split into the term n^hij and a contribution — un^ to the physical 
Hamiltonian density. Variation of the action with respect to h^ and tt %: * yields 

. f Till r rrff 

^ .(x) = - , = -4 . (50) 

n tor" (a) 5hij(x) y ' 

To evaluate these variational derivatives, one must use the relations 

Shij(x) ~ kl 3 . 57r lJ (x) 1 r i .-_fc7 S . 3/ , , ri s 

~ = <ffi<r(s - y) , = = — h %3 ii ki 5 A (x - y) , (51) 

and 

5 JhM - n - 8*&< x _ v ) ( 5 2) 

<^(y) _U ' 5^(y)- dkl6[X V) ' ^ 
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Defining the Poisson bracket appropriate to the reduced phase space, 

= r 3 / 5A 5B 5A 5B \ 

' J X \d%{xj 5¥i{x) ~ 6^(x)S~h ij (x) J ' 

any quantity A(hij,TT 13 ,t) obeys the equation of motion 

A=°— + {A,H"}. (54) 

Variation of the action with respect to N l yields the three constraints 

Hi ~ . (55) 

As before, the time evolution of N l is unconstrained by the action; in the absence of a 
gauge-fixing procedure, the three functions N % are arbitrary. 

3.4 Constraint Properties & Degrees of Freedom 

By lengthy direct calculation, it is possible to prove that the constraints Hi are first class, i.e., 
{Hi(x),Hj(y)} ~ 0. Furthermore, by applying the equations of motion to Hi, it is possible 
to show that Hi ~ 0, so the constraints are preserved by the equations of motion. We defer 
demonstrations of these two facts to section 5, where we will examine general relativity in 
the context of a class of realistic theories. This is an important consistency check, because 
a priori it is not clear that our procedure for solving the Hamiltonian constraint will yield a 
consistent action on the reduced phase space. 

As a final check, we revisit the counting of degrees of freedom in spatially covariant 
general relativity. After imposing cosmological gauge and solving the Hamiltonian constraint 
Ho ~ 0, general relativity is a theory of a unit-determinant spatial metric hij and its traceless 
conjugate momentum tt ij , so the theory contains ten canonical (or 5 real) variables. This 
reduction in the size of the phase space is compensated by a corresponding reduction in 
the number of constraints and arbitrary functions: the theory contains three first class 
constraints Hi ~ 0, and its equations of motion involve three arbitrary functions N % . 

5 • hij 's + 5-r j 's-3-Hi's-3- N if s = A canonical DoF . (56) 

Spatially covariant general relativity thus contains four canonical (or two real) degrees of 
freedom, the same number as fully covariant general relativity. 
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4 Ultralocal Modified Gravity 



We have two criteria in mind for our modified theories of gravity: two graviton degrees 
of freedom, and manifest spatial covariance. Our starting point is the action of spatially 
covariant general relativity, which has both of these properties. To modify general relativity, 
we will change the functional form of the scalar quantity 7r w , which in general relativity obeys 
= ttqr- This yields the action 

S = J dt d 3 x ( n^hij + W7r w - iV% ) , 
Hi = -2/^V fc 7F fc - uViit^ , (57) 

where is an unspecified scalar function of t, the phase space variables hij and n 1 ^ , and 
spatial derivatives. This action leads to the equation of motion 

8 A 

A = — + {A,H}, (58) 



where the Hamiltonian H is 

H 

and the Poisson bracket is 



J d 3 x ( -cun^ + N'H, ) , (59) 



{AB} = J d 3 x ^ 



5A 5B SA 5B 



Shij(x) Stt^(x) 5^{x)8h ij [ 



X 



(60) 



Retaining the manifest spatial covariance of the theory amounts to demanding 1) that the 
modified "Hj remain first class, i.e., 

{^(x),^(y)}~0, (61) 

and 2) that the modified constraints be preserved by the modified equations of motion, i.e., 

Hi~0. (62) 

Any theory satisfying these two points will be manifestly covariant under spatial diffeomor- 
phisms, with the constraints Irii acting as the generators of the gauge symmetry. Moreover, 
the presence of three first class constraints Irii on the phase space (hij,^) guarantees that 
such a theory contains two local degrees of freedom, exactly as desired. 
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In the remainder of the paper, we examine two classes of theories. First, for pedagogical 
purposes, we treat the case when is an ultralocal 7 function of time t and the phase space 
variables hij and 7p J ; in other words, 7r w will not contain spatial derivatives. Second, to 
make contact with general relativity, we treat the more realistic case when 7r w also depends 
on R, the Ricci scalar of h^. Forthcoming work will examine more general classes of scalar 
momenta [9] . In this section, we use the ultralocal case to introduce the formalism needed to 
determine when the constraints %i remain first class and when the constraints are preserved 
by the equations of motion. In section 5, we apply the formalism to the realistic case. In both 
the ultralocal and the realistic case, the consistency of the constraints with the equations 
of motion requires 7r w to satisfy an analogue of the renormalization group equation; scalar 
momenta satisfying this equation are manifestly invariant under rescaling of the volume 
factor oo. In the ultralocal case, this is the only consistency condition that arises. In the 
realistic case, demanding that the constraints %i satisfy a first class algebra is equivalent to 
demanding that 7r w obey a rather complicated differential equation. 

4.1 Constraint Algebra 

In this section, we will compute the Poisson bracket {'Hi{x),'H a (y)} assuming that 7r w is an 
ultralocal function, and use the result to determine when the constraints Hi remain first 
class. To simplify the calculation of {Hi(x), H-aiy)}, we split Irii into a tensor part Ji and a 
scalar part /Q. Concretely, we define the vector densities 

Ji = -2h l3 V k V k , Ki = -ujV^ , (63) 

in terms of which %i becomes simply 

H i = J i + lC i . (64) 

7 Local functions depend on the value of fields in a neighborhood around a point, so they can be written in 
terms of fields and derivatives of fields. Ultralocal functions only depend on the value of fields at a point, so 
they do not contain derivatives. The ultralocal limit is commonly used to analyze long distance cosmological 
perturbations. The idea of using the cosmological gauge in the ultralocal limit, sometimes referred to as the 
separate universes approach, is treated in [40] . 



18 



The Poisson bracket {Hi(x), / H a (y)} can then be written as the sum of more manageable 
brackets, 

{ni(x),H a (y)} = {Ji(x),J a (y)} + {Ki(x),K, a {y)} 

+ {J t (x), K a {y)} + {ld(x), J a (y)} . (65) 

To simplify the evaluation of these component Poisson brackets, we introduce the smoothing 
functionals 

Fj = J d 3 x fJi , F K = f d 3 x f% , 

Gj = J d 3 yg a J a , G K = J d 3 y g a IC a , (66) 

where the functions /' and g l are time-independent smoothing functions. We then compute 
the brackets 

{Fj,Gj} = J d 3 xd 3 yf i (x)g a (y){J i (x),J a (y)}, 
{Fj,G k } + {F k ,Gj} = J d 3 xd 3 yf i (x)g a (y)(^{J i (x) : }C a (y)} + {}C i (x) : J a (y)}y 

{Fk, Gk} = J d 3 xd 3 y f(x)g a (y){IC t (x),ICa(y)} ■ (67) 

We make the key assumption that the smoothing functions decay so rapidly at infinity that 
when we integrate by parts inside the smoothing functionals, the boundary term vanishes 
identically; the smoothing functions are otherwise arbitrary. With the freedom to integrate 
by parts at will, it is straightforward to compute variational derivatives of the smoothing 
functionals, and thereby to obtain explicit expressions for their Poisson brackets. By com- 
paring these explicit expressions to the formal expressions in equation (67), we will derive 
explicit expressions for the Poisson brackets involving Ji and /Q. 

To compute variational derivatives of the smoothing functional Fj, first integrate by 
parts to obtain 

Fj = 2 J d 3 xh ij 7t jk V k f\ (68) 

from which it follows that 

5Fj = J d 3 x [27r jk (W k f )5~h ij + 2~h ij (Vkf)5TT jk + 2~h ij n jk 5V k f i } . (69) 
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The first two terms in this integral are in a convenient form for taking variational derivatives 
with respect to and jt* k , but the third term requires finessing. To evaluate SV k f\ expand 
the covariant derivative as V k f l = d k f l + T' l kr f r ', where T l jk is the connection of the metric 
hij. It follows immediately that 5V k f % = f r $F kr . The identity 

61%. = h m (V r 5~h km + V k 8~h rm - V m 5h rk ) (70) 

thus implies that 2^ k h,ij5V k f l = f' i; K jk Vi5hj k , so equation (69) becomes 

SFj = J d 3 x {2^ k (W k f i )6h ij + 2h ij (y k f)6i? k + r^ k Vi6h jk } . (71) 
Integrating by parts, this reduces to 

5Fj = j d 3 x \2V k {V k r)5~h t] - Viifn^Shjk + 2h tJ (V k f l )5^ k ] . (72) 

From this expression, it is straightforward to compute variational derivatives of Fj, 

AF ~ 2 

' : ' ' = 2 ~5™ ^ k V k f - V, (f 7r mn ) - -n mn V l f , 



5ir mn 

The corresponding results for Gj are 



2€> y V,f. (73) 



Tr = 2( U T Vcfi- - V a (0 7T J — -7T V a Q , 



5n 



mil 



2 Cn KiN c g a ■ (74) 



The variational calculation for the smoothing functional F K is less straightforward. 
After integrating by parts, Fk becomes 



F K = to J d 3 x (dif) n u , (75) 

from which it follows that 

5F K = u j d 3 x (d t f) 8n u . (76) 

To evaluate 5tt w in full generality would be very difficult, so we will make some simplifying 
assumptions about the form of 7r w . In this section, we will assume that ir^ is an ultralocal 
function of t, h^, and n 1 *. 
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To facilitate calculations, we will enumerate all the scalars that can be built by con- 
tracting factors of against factors of . We begin by recursively defining H tj (n), the 
linked chain of n factors of 7t lj . The chain of zero factors of 7T lj is simply 

U ij (0) = h ij . (77) 

The process of adding a link to the chain is defined by 

rF(n+l) = r k U kj (n). (78) 

By closing the chain, one obtains scalars, 

0(n) = IT» . (79) 

The 4>(n) are the only scalars that can be built out of connected contractions of % and 7r tj . 
For an arbitrary ultralocal function ix w , it follows that 

°° Bit 

Since 0(0) = 3 and 0(1) = 0, 50(0) = 50(1) = 0. For n > 2, the variational derivatives of 
the 0(n) are 



5<f>(n)(x) _ ( ~ 5 rnn nIl{n) a b _ ^ mn n <j)(n - 1) ) 5 3 (x - y) , 



Sh mn (y) V 3 

~6^ n nU(n-l) ab 5 3 (x-y). (81) 



S(j)(n)(x) tab,,; . n 



5fr mn (y) 

The variational derivatives of F K are thus 
SF K 



CO 



5h 

rnn o 



W 5>^y \ 5 7^{ny k - f mn d>(n - 1)J , 



Similarly, the variational derivatives of are 



^ = w (a ( /<)y>^*n(»-%. (82) 



8h mn m=2 

m=2 ^ v 



21 



We emphasize that these results for F K and Gk rely on the ultralocality assumption, and 
will be modified in section 5. 

We are now in a position to compute the Poisson brackets of the smoothing functionals, 
from which we will extract the Poisson brackets of the vector densities Ji and /Q. 

• {Ji(x),J a (y)} 

To obtain the bracket {Ji(x) , J a {y)} , we first compute {Fj,Gj}. Combining the Fj 
and Gj variations into the bracket {Fj, Gj} yields 

{Fj, Gj} = zj d 3 z { (V c f ) (v^ a ) h ab n bc - (V k g a ) (V a f ) h l3 V k 

+ (V fe f ) V a (g a ~h t ^ k ) - (v c g a ) V, (f h ab n bc ) } . (84) 

After integrating by parts, using the definition Ji = —2hijVk& k , and using the identity 
(ViVj - VjVi) V a = R a bij V\ this reduces to 

{Fj, Gj} = J d 3 z {pj a V t g a - g a J$ a f + 2fg a V k [k jika + R jaik ) } . (85) 

From the symmetries of the Riemann tensor 8 and the traceless momentum tensor 9 , it 
follows that = V* (R ]lka + R jatk ), so the last term in the integrand vanishes. The 
connection terms inside the remaining covariant derivatives cancel to yield 

{Fj, Gj} = J d?z (f l J a d ig a - g a Jid a f) . (86) 

To extract the bracket {Ji{x) , J a (y)} from this result, first relabel dummy indices 

{F J? Gj} = j d?x fJad l g a - J d 3 y g a Jid a f . (87) 

Under the spatial derivatives in this equation, insert the identities 

g a (x) = J d 3 y 6 3 (x - y)g\y) , f (y) = J d 3 x 6 3 (x - y)f\x) , (88) 

to obtain 

{Fj,Gj} = J d z xd?yf{x)g a {y) (j a (x)d xi 5 3 (x - y) - J^d^x - y)) . (89) 



Rabcd — Rcdab-i Rabcd Rbacd Rabdc- 
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Comparing this expression to equation (67) yields the identity 

{Ji(x), J a (y)} = J a (x)d xi 5 3 (x -y)- J % {y)d y «5\x - y) . (90) 

This is the same algebra obeyed by the Hi in equation (15). This result is completely 
independent of our choice of n^, and will carry over unchanged into section 5. 

To obtain {Ji(x),K, a (y)} + {/Q(x), J a (y)}, we first compute {Fj,G K } + {F K ,Gj}. 
Assembling the Fj and Gk variations into the Poisson bracket {Fj, Gk} yields 

/°° f) 
d 3 z (d a g a ) d^t) mYi{m ~ ■ ( 91 ) 

By expanding the covariant derivative, simplifying the ensuing total derivative of 7r w , 
and recalling that /Q = — wVi7r w , this expression reduces to 

d 3 ^ r IC t d a g a -uj Jd 3 z (dif) (d a g a ) d^n) m ^ m ) ■ ( 92 ) 

Similarly, 

d 3 z g'ICadif + u Jd 3 z (dif) (d a g a ) £ ^m^(m) , (93) 

so the sum of the two brackets simplifies considerably, 

{Fj, G k } + {F K , Gj} = J d 3 z (f}C l d a g a - g a JC a d t f) . (94) 
Integrating by parts and invoking the identity <9j/C a = <9 a /Q yields 

{Fj, G k } + {F K , Gj} = J d 3 z (fJC a d t g a - g^dj 1 ) . (95) 

To extract the quantity {^(x), K, a (y)} + {/Q(:r), J a {y)}, relabel dummy indices and 
insert the identities in equation (88) to obtain 

{Fj,G k } + {F k ,Gj} = J d 3 xd 3 yf(x)g a (y) 

x {K, a {x)d xi 5\x -y)- Ki{y)d^{x - y)) . (96) 

Combined with equation (67), this result implies that 

{Ji(x), K a {y)} + {Kiix), J a (y)} = K a {x)d xi 8 3 {x -y)- K i {y)d ya 8 3 {x - y) . (97) 

This expression depends strongly on the assumed form for 7r w . This result is modified 
heavily in section 5.1, when is allowed to depend on R. 
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• {/Q(x),/C a (y)} 

To obtain {/Q(x), /C a (?/)}, we first compute the bracket {F k ,Gk}- Substituting the 
F K and Gk variations into the Poisson bracket {F K , Gk} yields 

OO OO r, 

{F K ,G K } = u, 2 (d t f) (d a g a ) fw\ 

t^it^i d<P{m) d(j){n) 

x (Tl(n) bc Tl(m - l) bc - U( m y k U(n - l) jk ) . (98) 

From the definition of the momentum chain H{n)^ , it follows that U(n) bc U(m — 1)^ = 
Yi{my k Ii{n — = <f>(n + m — 1). The terms of the sum thus vanish order by order, 
so the bracket reduces to 

{F K ,G K } = 0. (99) 

By comparing this result to equation (67), it is apparent that 

{JCi(x),K a (y)} = 0. (100) 

When 7r w is an ultralocal function of the phase space variables, the Poisson bracket 
{K,i(x), )C a (y)} vanishes identically. This will not be the case when depends non- 
trivially on R, as in section 5.1. 

By substituting equations (90), (97), and (100) into equation (65), and recalling that T-Li = 
Ji + K,i, we obtain 

{Hi(x),Hi(y)} = Hj(x)d xi 5 3 (x -y)- U t (y)d y >5 3 (x - y) . (101) 

This is the same algebra obeyed by the Hi in equation (15), and by the Ji in equation (90). 
Since Hi ~ 0, this result implies that {Hi(x),T-Lj(y)} ~ 0, so the constraints Hi are first 
class. To establish this result, we assumed only that 7r w was an arbitrary ultralocal function 
of t, hij, and fc^; we showed that this was equivalent to making 7r w a function of t and the 
scalars (f>(n) defined in equation (79). Evidently, 7r w can be made any ultralocal function of 
the phase space variables and the momentum constraints will remain first class. 



4.2 Consistency of Constraints with Equations of Motion 

In this section, we will compute the time derivative Irii assuming that 7r w is an ultralocal 
function, and use the result to determine when the constraints Hi are preserved by the 
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equations of motion. The time evolution of Hi is determined by the equation of motion 

■ dH 

Hi = -^ + {Hi,H}, (102) 



where 

H 



J d 3 x ( -un u + WiU ) . (103) 



Since Hi = Ji + /Q and dJi/dt = 0, it follows that &Hi/dt = dICi/dt. Recalling that 
Ki = —udiTTu, the first term in equation (102) becomes 

^ -a V (io4) 



0i ' V 9t , 

To simplify the bracket {Hi,H}, we define 

IT, = J d 3 ^, (105) 

so that iJ can be written as 

H = -uIL,, + / d 3 x WHi . (106) 



From the first class character of the constraints Hi, it follows that {Hi,H} ~ —u{Hi,H u }. 
Since Hi — Ji + JCi, the second term in equation (102) becomes 

{Hi, h} ~ -w{ ji, ru - w{/c 4 , ru . (107) 

To compute the brackets {J7i, U^} and {/Q, IT,}, we first compute the smoothing functional 
brackets 

{F J ,U U } = J d^fWI^fx),^} 

{F K ,U U } = J d 3 xf(^){^(x),n w }. (108) 

We have already done all the work needed to evaluate these two brackets: since Il w can be 
obtained from Gk by the substitution d a g a — > brackets involving n w can be obtained 
by applying this substitution to brackets involving Gk- 
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{JiM 

To compute the bracket {^,11^}, we first compute the bracket {F J: n w }. Applying 
^a5 ,a — > ^ 1 to equation (92) and integrating by parts yields 

{Fj,^} = j d 3 xfd t + f^m^m)-^-^ . (109) 

It follows by comparing this result with equation (108) that 

/ oo d \ 



m=2 



• {/Ci,n w } 

To compute the bracket {/Q, 11^}, we first compute the bracket {F K , U u }. By applying 
the transformation d a g a — > lu^ 1 , equation (99) becomes 

{F K ,U 1tl } = 0. (Ill) 

Along with equation (108), this implies that 

{/Q,ru = o. (112) 

By substituting equations (110) and (112) into equation (107), we obtain 

Upon inserting equations (113) and (104) into the equation of motion (102), the ud^^ terms 
cancel to yield 

*~-*( u £ + *?™« m >fl&)- (114) 



m=2 

Demanding % ~ implies the consistency condition 

<-\ OO <-\ 

07T„, . ^ j / \ 07T (I 
m=2 

where f(t) is an arbitrary function of time. We observe that the equation of motion (58) is 
invariant under 7r w — > n u + g(t), where g(t) is an arbitrary function of time, so we are free 
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to apply this transformation to simplify our consistency condition. If we choose g{t) so that 
cog'(t) = f(t), the consistency condition becomes 

dir ()';% 

u W + *^^ m) dfifi~°- (116) 

m=2 rv 1 

By assumption, uj(t) is an invertible function of time, so d/dt = uid/dui. Our consistency 
condition can thus be written as 

Avr^-O, (117) 

where we have defined the operator 

d <) 

A ^ + 5 2 m0(m to' (118) 

To rule out the possibility of a 7r w which satisfies An^ ~ while Air^ ^ 0, we note that 
the constraints %i contain one power of spatial derivatives, while by assumption the scalar 
momentum is ultralocal. To satisfy An u ~ 0, the quantity An^ would need to depend on 
the constraints and would thus need to contain at least one power of spatial derivatives. 
However, applying A to ir^ does not increase the number of spatial derivatives. It follows that 
Anu cannot contain any spatial derivatives, and thus cannot depend on "Hj. The consistency 
condition can therefore be promoted to 

Att w = 0. (119) 

To obtain the most general solution to this equation, we first note that A (w _n 0(n)) = 0, 
which motivates us to define 

0(n) = *™ . (120) 

The most general solution to the condition A-n^ = is an arbitrary function of the 4>(n). 
The explicit time dependence of is thus determined by its dependence on the phase space 
variables. 

To understand this result, we return briefly to the phase space (hij,ir l: >). To construct 
three-scalars out of the tensor and the traceless tensor nji, we begin by recursively 
defining n^(ra), a chain of n factors of 7f^ linked together by factors of h^. In analogy with 
our construction of the 4>{n) of equation (79), we define 

ITi'(O) = h ij = Q-^iO) , (121) 
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and 

Il?(n + 1) = r T a h ab U b \n) = u-^UlFin) , (122) 

from which it follows that U l ^(n) = VL~ l uo^Wi (n) . The contraction hijlCp(n) yields the 
desired scalars, 

Mn) = ^{n) = ^ . (123) 

The <j>r(n) are the only scalars that can be built out of fully connected contractions of hij 
and 7Ty. In the presence of the constraint ou ~ w(t), it follows that 

T (n) ~ 0(n) . (124) 

In other words, the 4>(n) are the scalars on the phase space (h^, n lj ) which have the correct 
conformal weight to have been derived from three-scalars on the phase space (hij,ir l: >). It 
follows that the 0(n) are invariant under a rescaling oj — > fiu of the volume factor u>, and 
the condition An^ = is thus analogous to a renormalization group equation. 



4.3 Summary 

In this section, we developed a formalism for testing when our modified theories of gravity 
lead to a consistent first class constraint algebra, and hence contain two degrees of freedom. 
To develop the formalism, we made the simplifying assumption that the scalar momentum 
7r w is an ultralocal function of time t and the phase space variables and 7p J . This as- 
sumption is sufficient to guarantee that the constraints T-Li remain first class. However, for 
the constraints to be consistent with the equations of motion, 7r w must be invariant under 
renormalization of the volume factor u. Concretely, 7r w must obey the renormalization group 
equation 

Att w = 0, (125) 

where 

d <) 
A^- + 5># m) — . (126) 

m=2 rv ' 

Satisfying this equation completely fixes the dependence of ir u on u>(t). In the next section, 
we will apply the methods of this section to generalize this result to a more realistic class of 
scalar momenta. 
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5 Realistic Modified Gravity 



The ultralocal ansatz has the virtue of simplifying calculations, but it has the defect of being 
manifestly unphysical: the laws of nature are local, not ultralocal. In this section, we will 
apply the formalism developed in the last section to theories in which 7r w depends on spatial 
derivatives of the metric hij through a dependence on the Ricci scalar R. Since the 7Tgr of 
spatially covariant general relativity belongs to this class, we call it the "realistic" class. As 
we will demonstrate, realistic 7r w must obey stringent consistency conditions in order for the 
Hi to generate a consistent first class constraint algebra. 

5.1 Constraint Algebra 

In this section, we will compute {'Hi(x),'H a (y)} assuming that ir^ is a function of t, the phase 
space variables and tt^, and the Ricci scalar R. We will then use the result to determine 
when the constraints Irii remain first class. 

As before, we decompose "H, into a tensor part Ji = —2hijV and a scalar part 
/Cj = — wVi7r w . Computing {Hi(x), % a (y)} is then a matter of computing the four brackets 
in equation (65). The result for {J7i(x), J a {y)} carries over unchanged from equation (90), but 
we will have to revisit the brackets involving /Q. To do so, we will first evaluate the smoothing 
functional brackets {Fj, G K } + {F K , Gj} and {F K , G K }. By comparing the ensuing explicit 
expressions to the formal expressions in equation (67), we will derive explicit expressions for 
the Poisson brackets involving /Q. 

Our analysis of the variational derivatives of the smoothing functional Fk defined in 
equation (66) proceeds exactly as in the ultralocal case up to equation (76), where the 
quantity 5-k w arises. In this section, we assume that 7r w is a function oft, h^, , and R. To 
simplify calculations, note that this is equivalent to making 7r w a function of t, R, and the 
<f>(n) defined in equation (79). It follows from this assumption that 




(127) 



Substituting this result into equation (76), using the identity 5R = —R? k 5hjk + V k V^5h 
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and integrating by parts yields 

5F K = uJ d 3 , ( di f) (± -jjjfcsm - |p^ N 

+u J d 3 * V^V fc (jdif) ||) 6h jk . (128) 
Using equation (81), it is now straightforward to compute the variational derivatives of F K , 

-co {dtf) ^~5JTR Jk + co~SJT^ k (W) ||) , 
i^L = u ( d .f) y nj^sg nfn _ i) . (129) 

The corresponding results for Gk are 

. „ „. v d'Ku, ( ~„„, , 1 



3h mn , 



pi / 1 N 

E (C n n(m) bc - -^ m >(m - 1)^ 

-a; (0^) ^5lTR bc + u~5™V b V c (W a ) ||) , 

= uj (d a g a ) y m-^-^ c n n(m - l) bc . (130) 
<57T mn ^ <90(m) mn v 7 v ' 

We are now in a position to compute the brackets involving /Q. 

• { t 7i(x),/C a ( 2 /)} + {/C i (x), t 7 a ( S /)} 

To compute { ^(x), /C a (y)} + {/Q(z), J a (y)}, we first compute {Fj,G K } + {F^Gj}. 
We begin by substituting equations (73) and (130) into the bracket {Fj,G K }. After 
expanding and simplifying a total derivatives of (f>(n), {Fj,Gk} turns into 



{Fj, G K } = -u> f d*z P (d a g a ) ^rV^(m) + 2oj [ d 3 z (d a g a ) ^R^kf 
-2u J d 3 z (V fe f) V,V fe ((d a g a ) C - 
+| W J d s z (dif) V C V C ((d a g a ) |j 



/ d -z {dif) W) ^ + |>M^J • (13D 
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To finesse this expression, integrate by parts, use the identities VjV 'jV l = V jViV 1 + 
RijV 1 and 2V jR^ = Vj-R, simplify a total derivative of n w , use the identity /Q = 
— cuViiTu, and expand to obtain 

chr u 



{Fj,G K } = J dhfKidag* + d 3 zV k (dif) V fc {{d a g a ) 



dR 



-u, / d», (a,/-) W ) (§*f| + E-^H^R ) • < 132 > 

Similarly, 



+c / d 3 , (ft/*) (0^) fef| + f>*(m)^) , (133) 



so the sum of the two brackets reduces to 



{Fj,G k } + {F k ,Gj}= J d 3 zf i K i d a g a - J d 3 zg a K a dif 
+^co I d 3 zV k (dif) V fe (W)f) 



~u J d 3 zV k (d a g a )V k (faf) ||) . (134) 
After integrating by parts, expanding, and using the identity diK a = <9 a /Q, this becomes 
{Fj,G k } + {F k ,Gj} = J d 3 zf]C a d i g a - J d 3 zg a K i d a f 

+ J d 3 z {dif) (d k d a g a ) M k - J d 3 z (d a g a ) (dAf) M k , (135) 

where 

^* = -* wV *7§- ( 136 ) 
o OR 

To extract the bracket { Ji(x), JC a (y)} + {lCi(x), J a {y)}-> integrate by parts, relabel 
dummy indices, and insert the identities in equation (88) to yield 

{Fj, G K } + {F K , Gj} = J d 3 x d 3 y f(x)g a (y) {K a {x)d xi 5 3 {x - y) - K i {y)d y a5 3 {x - y)) 

+ / d 3 xd 3 yf(x)g a (y)d xi (-M k (x)d xk d x aS 3 (x - y)) 

d 3 xd 3 yf{x)g a {y)d r (-M k (y)d yk d yi S 3 (x - y)) . (137) 
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By comparing this expression to equation (67), it is clear that 

{Ji(x), K a {y)} + {Ki{x), J a (y)} =)C a (x)d xi 5 3 (x - y) - Ki(y)d y «5\x - y) 

+d x i {—M k (x)d x kd x a8 3 (x - y)) 
-d ya {-M\y)d yk d t 5\x - y)) . (138) 

• {Ki(x),K a (y)} 

To compute {fCi(x), fC a (y)}, we first compute the bracket {F k ,Gk}- Substituting 
equations (129) and (130) into the bracket {Fk, Gk} yields 



{F K ,G K }=u 2 J d 3 z (d a g a ) ^V^V^(^r) 



OR 



J d '*M^ s * k {(W)%%), d39) 



where 

n=2 

After integrating by parts and expanding, the bracket becomes 



|^ = $E»^n(„-i) fc . (140) 



{Fk, G k } = J d 3 z (d,f) (BAg") W - j i 3 z (,%g") (d k d t P) N k , (141) 

where 

To extract the bracket {K,i(x), K. a (y)}, integrate by parts, relabel dummy indices, and 
insert the identities in equation (88) to obtain 

{F K ,G K }= f d 3 xd 3 yf(x)g a (y)d xi (-M k (x)d xk d x a5 3 (x - y)) 

-J d 3 xd 3 yf i (x)g a (y)d ya (-Af k (y)d yk d yi 5 3 (x - y)) . (143) 

Comparing this expression to equation (67), it follows that 

{1d(x), K a (y)} =d xi (-N k {x)d xk d x a5 3 (x - y)) 

-d ya (-Af k (y)d yk d yi 5 3 (x - y)) . (144) 
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By substituting equations (90), (138), and (144) into equation (65), and recalling that Hi = 
Ji + /Q, we obtain the identity 

{Uiix^iijiy)} =H j (x)d xi 5 3 (x-y)-H i (y)d yj 8 3 (x-y) 

+d xi (-X k (x)d xk d xj 8 3 (x - y)) - d yj (-I h (y)d yk d yi 8 3 (x - y)) , (145) 

where X k = Mk + A4, or 

Si? Ott^* dvr^ Si? 3 Si? 
Expanding the derivatives in this expression and using the fact that Irii ~ yields 

{H t {x),H 3 {y)} ~ - (Z fe (a;) + X fe (y)) d^d^x - y) 
-{d xi I k {x))d xi d xk 5\x-y) 

+ (d yj X k (y))d yi d yk 5 3 (x-y). (147) 

The three terms of this equation are algebraically independent, so the necessary and sufficient 
condition for the Poisson bracket {i-Li(x),l-Lj(y)} to vanish is 

X fe ~0. (148) 

In the ultralocal case the constraints were automatically first class, but to generate a first 
class constraint algebra in the realistic case, the scalar momentum must obey the fearsome 
looking differential equation Xk ~ 0. 

As a check, we will now compute the Xk arising from the 7Tgr of spatially covariant 
general relativity. Recall from equation (46) that 



n GR = -W-y'w- 2 0(2)-^- 2 / 3 i? + 2A. (149) 
Since 7r GR is a function only of 0(2) and R, its partial derivative with respect to simplifies, 

-d^ = 2 dm 7Tjk - (150) 

After substituting this relation into the definition of Xk and recalling that Ji = — 2hij'Vk^ k , 
Xk becomes 

. 2 ~ ( d^GR ™ <9vr GR <9tt gr ~ , dir GR \ 
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Upon substituting the derivatives 

dn GR 4 1 <9tt gr 



90(2) 3w 2 7r GR ' OR " 3W 2 /3 7T GR 



(152) 



into I/c(7Tgr), the term in parentheses vanishes. By using the relations /Q = — wVj7r w and 
Hi = Ji + ICi, we obtain 

16 

Xfc( " GR) = 9^7%|^^- ( 153 ) 

Since Hi ~ 0, the scalar momentum 7Tgr satisfies ~ 0. The constraints "Hj of spatially 
covariant general relativity thus generate a first class constraint algebra. 



5.2 Consistency of Constraints with Equations of Motion 

In this section, we will compute the time derivative Hi for realistic 7r w assuming that the 
constraints Hi are first class, and use the result to determine when the constraints Hi are 
also preserved by the equations of motion. The analysis of Hi proceeds exactly as in the 
ultralocal case until we arrive at the expression 

4 = ~di (uir u + u^j - ou{Ji, IU - cu{}C t , IU , (154) 

where as before 




(155) 



The point of departure from the ultralocal case is the evaluation of the two Poisson brack- 
ets {jZjjIXj} and {/Cj,]!^}. To compute them, we first compute the smoothing functional 
brackets{.Fj, 11^} and {Fk, IT,}. As in the ultralocal case, we will obtain brackets involving 
by applying the substitution d a g a — > u^ 1 to brackets involving Gk- 

• {JiM 

To obtain the bracket {J7i,n w }, we first compute the bracket {Fj,!!^}. Applying 
d a g a — > to equation (132) and integrating by parts yields 

{Fj, IU = J d 3 x f% (-k u + 2 -R^ + f; W>M^y - ^V k M k j , (156) 
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where as before 

4 ~ 8-k 

M k = --uV k -^. (157) 
It follows from an application of equation (108) that 

(Of) 00 f) \ 

-t w + -=R^g + V m<f>(m)— ^- - u;- 1 V fc M fe . (158) 

• {/Ci,n w } 

To compute the bracket {/Q,]!^}, we first compute the bracket {F^n^}. After sub- 
stituting d a g a — >■ and integrating by parts, equation (141) becomes 

{F K , IL} = J d 3 x f% {-u- 1 V k U k ) , (159) 



where as before 



JV4 = W ^ v 9^- w »ff v aF (160) 

Comparing with equation (108) yields 

{K i ,n u } = d i (-u>- 1 v h N' k ) . (i6i) 

After substituting equations (158) and (161) into the equation of motion (154) and recalling 
that X k = M. k + A/fc, we obtain 



4 = -ft J^- + + W>(™)£^r " *<" _1 V*Z* • (162) 

V at 3 Si? ~ <90(m) / 



m=2 

Since the Hi are assumed to be first class, it follows necessarily that X k ~ 0. Demanding 
Hi ~ thus implies the consistency condition 

c>t 3 OR ~ 2 a0(m' 



where fit) is an arbitrary function of time. Recall once again that the equation of motion (58) 
is invariant under n u — > ix w + g(t), where g(t) is an arbitrary function of time. By choosing 
a function g(t) such that oug'(t) = fit), the consistency condition becomes 

w— V -wR — — + w> mSim)-^—. — -~0. (164) 
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Since u(t) is assumed to be an invertible function of time, d/dt = ud/du. In analogy with 
our approach in the ultralocal case, we rewrite the consistency condition as 



Att w ~ , (165) 
where we have redefined the operator A as 



AE/ + -^ + y m ^(m)-f-. (166) 
du 3 OR ^ 9 1 >d^(m) 1 ; 



To rule out the possibility of a ir^ which satisfies Att^ ~ while A7r w ^ 0, we note that the 
constraints Hi contain a term Vi7r w , making the constraints higher order in spatial derivatives 
than 7r w itself. However, by examining a series expansion of 7r w in the parameter R, one can 
verify that applying A to does not alter its order in spatial derivatives. 10 It follows that 
AtToj cannot depend on Irii. The condition Air u ~ is therefore equivalent to the apparently 
stronger condition 

Att w = 0. (167) 
Since A(w~ n 0(n)) = and A(u~ 2 / 3 R) = 0, we are led to define the quantities 

The most general solution to the condition Att^ = is an arbitrary function of R and the 
4>(n). In this manner, the dependence of 7r w on oj(t) is determined by its dependence on the 
phase space variables. 

As before, to understand this result, we return briefly to the phase space (hy, 7r y ). As 
shown in section 4.2, the only scalars that can be built out of the tensor hy and the traceless 
tensor ttj, are the 4>t(h) — oo~ n (j)(n). If we impose the gauge-fixing constraint uj ~ ou(t), then 
4>T{n) ~ (f>(n); likewise, the Ricci scalar R of the metric h^ obeys R ~ R. 11 This means that 
R and the 4>(n) have the correct conformal weight to have been derived from three-scalars 
on the phase space (hij, n 1 ^). The scalars R and the <t>(n) are thus invariant under a rescaling 



10 Spatial derivatives enter 7r w solely through R, so the derivative expansion of can be written 7r w = 
Y^k=Q c kR k i where the coefficients depend on co and the <p( n )- Applying the A operator to 7r w changes 
the functional form of the Cfe, but does not generate higher order powers of R. 

11 See equation (259) in appendix C. 
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uj — > fiuj of the volume factor uj, so once again An^ = is revealed to be analogous to a 
renormalization group equation. 

As a check, we will now apply the renormalization group equation to the scalar mo- 
mentum 7Tgr of spatially covariant general relativity. Since 

7rGR = _ vf \M 2 ) _i?+2A ' ( 169 ) 

the scalar momentum ttgr satisfies the condition A^gr = 0; this implies that the constraints 
of the theory are preserved by the equations of motion. Combined with the result that 
X k (7iGR) ~ 0, which implies that the constraints are also first class, it is now clear within 
the context of our formalism that the constraints %i of spatially covariant general relativity 
generate a consistent first class algebra. This result justifies the assertions we made in the 
first paragraph of section 3.4. 



5.3 Summary 

In this section, we applied the formalism developed in section 4 to determine when scalar 
momenta 7r w built out of hij, Tr^, and R yield a consistent first class constraint algebra. To 
ensure the first class character of the constraints Hi, it is necessary and sufficient for n u to 
obey the condition 

Z fe ~0, (170) 

where 

Xk = u ^7f7 V TFlk ~ u TFlk V ^7f7 ~ o uV k—~- ■ 171 
OR 07T J * on^ OR 3 dR 

If dix^jdR = 0, then X k = 0, so ultralocal scalar momenta satisfy this condition trivially. 
The scalar momentum 7Tgr of spatially covariant general relativity depends essentially on R, 
and thus satisfies this condition non-trivially. 

To guarantee the preservation of the constraints Irii by the equations of motion, the 
scalar momentum tt^ must also be invariant under renormalization of the volume factor uj. 
This requires 7r w to obey the renormalization group equation 

Att w = 0, (172) 
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where 

o Of) °° f) 

A = u-^- + + Y m<P(m)-£— . (173) 

duo 3 ^ ; <90(m) 

This is a generalization of the renormalization group equation (125) to include a possible 
dependence of 7r w on R. The scalar momentum 7Tgr satisfies this condition in addition to 
the first, so the constraints of spatially covariant general relativity generate a consistent first 
class constraint algebra. 



6 Conclusions 

In this paper, we developed a general formalism for verifying the consistency of spatially 
covariant modified theories of the transverse, traceless graviton degrees of freedom. It was a 
long road, so it is worth retracing our steps to see the logic of our path. 

In section 2, we showed how to express general relativity as a theory of a spatial metric 
hij and its conjugate momentum 7r y . In this language, the general covariance of the theory 
is represented on the phase space (/ijj, 7r u ) by the algebra of the four constraints In 
section 3, we showed how to obtain a spatially covariant version of general relativity. We 
began in section 3.1 by splitting the phase space (hij, ir^) into the phase space (ui, tv u ) of the 
spatial volume factor and the phase space (hij,^) of the unit-determinant metric. In the 
context of cosmology on an FRW background, it is natural to represent time diffeomorphism 
symmetry on the phase space (u, tt^) and to represent spatial diffeomorphisms on the phase 
space (hij,^); in section 3.2, we showed how to achieve this splitting using a cosmological 
gauge condition. On an expanding background, oj drops out of the dynamical phase space of 
the theory, and its conjugate momentum n u becomes the scalar Hamiltonian density on the 
phase space (hij,^); in section 3.3, we showed how to reduce the phase space by solving 
the Hamiltonian constraint in cosmological gauge. By successfully projecting the degrees of 
freedom of general relativity onto the reduced phase space (hij, 7f u ), we have shown how to 
represent the graviton dynamics of general relativity on the class of conformally equivalent 
spatial metrics. 

To modify general relativity, we simply modified the functional form of the scalar 
momentum 7r u while retaining the explicit spatial diffeomorphism symmetry generated by 
the three constraints "Hj. In section 4, we considered the case in which is an ultralocal 
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function of the phase space quantities hij and rc % K In this case, the consistency of the 
constraints Irii imposes a single non-trivial condition on the form of 7r w , namely that it 
must satisfy a renormalization group equation with flow parameter oj. The renormalization 
group equation encodes the fact that 7r w must be invariant under flow through the space of 
conformally equivalent spatial metrics. In section 5, we applied our formalism to the more 
realistic case in which is also allowed to depend on R, the Ricci scalar of the metric 
h^. In this case, 7r w must satisfy a corresponding renormalization group equation, but its 
form is further constrained by a differential equation that relates its dependence on R to its 
dependence on the phase space variables and 7r u . 

As a proof of principle, this paper demonstrates the possibility of consistently modifying 
the graviton equations of motion, but more remains to be done. In forthcoming work [9], we 
will apply our formalism to search for realistic alternatives to general relativity. In particular, 
we will examine scalar momenta with a more general dependence on derivative quantities, 
such as the Ricci tensor Rij of the unit-determinant metric h^, which fully determines the 
spatial curvature in cosmological gauge. After deriving consistency conditions in the more 
general case, we will attempt to solve them perturbatively to obtain valid scalar momenta 
7r w related to the 7Tgr of general relativity by the deformation of a continuous parameter. If 
we discover non-trivial modifications of general relativity that contain only two degrees of 
freedom, it could open up new lines of theoretical and experimental research. A null result, 
on the other hand, would serve as further evidence of the uniqueness of general relativity. It 
will be interesting to see just how far we can push this program. 
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Appendix A: Covariant Constraint Algebra of GR 



Recall the Poisson bracket of GR 

[ 3 / SA SB 5A SB \ 

1 ' 1 = J \Sh mn (z) ^ mn {z) ~ 6n™(z) 5h mn (z)) (7) 

and the constraints 

n = -VK(r - 2A) + -L _ 

Hi = -2h ij V k n ik . (175) 

Our object in this section is to derive the constraint algebra 

{n (x),n (y)} = W(x)d xi S 3 (x -y)- W{y)d yi S\x - y) 
{n (x),n i (y)}^n (y)d xi S 3 (x-y) 

{ni(x),nj( y )} = n,(x)d x ,s 3 (x - y )- n i ( y )d yi s 3 {x - y ) . (176) 

To evaluate these Poisson brackets, we first define the smoothing functionals 
F H = J d 3 xf(x)n (x) , F = J d 3 xf(x)U t (x) , 
G H = J d 3 yg°(y)Ho(y) , G = J d 3 y g\y)U a (y) , (177) 

where the functions /°, /*, g°, and g l are time- independent smoothing functions. We then 
compute the brackets 

{F H ,G H } = J d 3 xd 3 yf(x)g (y){Ho(x),Ho(y)} 

{F H ,G} = J d 3 xd 3 yf(x)g a (y){H (x),n a (y)} 

{F,G} = J d 3 xd 3 yf i (x)g a (y){n i (x),n a (y)}. (178) 

As in section 4.1, we assume that the smoothing functions decay so rapidly that they elim- 
inate all boundary terms generated by integration by parts, but that they are otherwise 
arbitrary. This greatly simplifies the explicit evaluation of the brackets of the smoothing 
functionals. By comparing the explicit forms of the brackets to the implicit forms in equa- 
tion (178), we will derive explicit formulae for the brackets of the T-L^s. 

40 



To simplify the calculation of the variational derivatives of F H , we will split the Hamil- 
tonian constraint "H m to a kinetic piece Ht and a potential piece Hy. Explicitly, we have 
Hq = Ht + "Hy, where 



H v = -Vh(R - 2A) . 



kl I ^7C kl 



Similarly, Fh = Ft + Fy, where 



F T = / d 6 xf(x)H T (x) , F v = d d x f(x)H v (x) 



Computing the variation 5F T is straightforward: 

5F T = J d 3 xf (1= (2tt>^ - 7r* fc 7r«) - \n T h^j Sh z 



It follows that 



5F T 

5Ft_ 

5n mn 



( 1 1 

= J \-y= ( 27r k^ k K ) - -n T n 



f ^yj- (27r mri h mn Ti ^ . 



Likewise, 



5Gt 

5Gt 
5ir mn 



Vh 



(2tt 



u7T — 7T .7T 



fl 1 (27T mn h rnn 7l fe ) 



;i79) 



;iso) 



(181) 



(182) 



(183) 



Keeping in mind that <5i? = —5hijR 1 ^ + V^V % Shij — V k V kh % ^ 5hij, computing <5Fy is 



just as straightforward: 



+ y d 3 x v^/° (V k V k h ij 5h tJ - V'VShij) 



(184) 
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Before taking variational derivatives, we exploit our freedom to integrate by parts to pull 



the covariant derivatives off the metric variation 5hif 



It follows that 

5F V 



5F V 



5ir mn 

Similarly, 

5G V 



SF V = J d 3 xf° (]^ivh lj + VhR ij ^j 5h i3 

+ J d 3 xVh V fc V fc /° - W/°) Shij . (185) 

/° (^K v h mn + VhR mn ^j + Vh (h mn V k V k f - V m V n f) 
0. (186) 



g° (\u v h mn + VhR mn J + \fh (h mn V k V k g° - V m V n g°) 
oh mn \2 J 

|^ = 0. (187) 

Before computing 5F, we integrate by parts inside F: 

F = 2 J d\rh,;^T,f' . (188) 
This simplifies the variational calculation: 

5F = 2 J d 3 x(y k f i )7r jk 5h lj + 2 J d 3 x(V k f )h ij 5^ k + 2 J d 3 x7r jk h lJ 5V k f l . (189) 

To evaluate 5V k f \ first expand the covariant derivative as V k f % = d k f l + T l ka f a . It follows 
that SV k f = f a 5T{ a . The identity 

Sr l ki = h lm (V t 5h km + V k 5h im - V m 5h lk ) (190) 

implies that 2ir jk hij5V k f l = f % ^ k V i5hj k . Substituting this result into the expression for 5F 
and integrating by parts yields 

5F = 2 J d 3 x iy k f)^ k 5h l3 - J d 3 x Vi(r ir jk )5h jk + 2 J d 3 x (Y,.f )h,/^ k . (191) 



It follows that 

5F 



Sh, 



2(V k f)n^ n -V l (fV 



SF 

^ =2{V k r)h^ k n . (192) 



5li mn 
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Likewise, 



5h 



<)G =2(V c s a )/wtV (193) 



5n mn 

We are now in a position to compute the Poisson brackets of interest. 

To calculate Hj(y)}, we will first calculate {F, G}. Substituting equations (192) 

and (193) into the Poisson bracket yields 

{F, G} = 2 J d 3 z h ab n bc (V c f ) {V l9 a ) - 2 J d 3 z ^-tH* (V^ a ) (V f ) 

- 2 y d 3 ^ {V c g a ) V, (r^Tr 60 ) + 2 J d 3 z (V fc r) V a (<7<%V fe ) . (194) 

After integrating by parts, applying the identity (VjVj — Vj Vj) -u a = R a bij u b , and 
recalling that "Hj = — 2/i ij -Vfc7r- 7 ' fc , this bracket becomes 

{F,G} = y d 3 ^ {fU a V ig a -g a UNaf l ) 

+ 2 y d 3 z V fc (% fea + i? jaife ) . (195) 

It follows from the symmetry (-R afecrf = Rcdab) and antisymmetry (i? a6cc z = -Rbacd = 
—Rabdc) properties of the Riemann tensor that Rj a ik = —Rkija- The symmetry property 
(7r u = 7T jl ) of the momentum tensor then implies that ir jk (Rjika + Rjaik) = 0, so 

{F, G} = y d 3 z {rU a V l9 a - g a U t V a r) . (196) 

Upon expanding the covariant derivatives, the connection terms cancel, yielding 

{F, G} = y d 3 z (m a d t g a - g a Uid a f) . (197) 

To extract the Poisson brackets {'H i (x),'Hj(y)}, first relabel integration variables, 

{F,G} = J d 3 x f(x)n a (x)d xi g a (x) - y d^g^n^dyafiy) , (198) 

then use the identities 

g a {x) = J d 3 y 5 3 (x - y)g a (y) , f{y) = J d 3 x 5 3 (x - y)f{x) , (199) 
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to write 

{F,G} = J d 3 xd 3 yf(x)g a (y)(n a (x)d xi S 3 (x-y)-n i (y)d y a5 3 (x-y)). (200) 
By comparing this expression to (178), we obtain the identity 

{Ui{x),nj{y)} = nj(x)d xi 5 3 (x -y)- Hi(y)d yj 8 3 (x - y) . (201) 
{Ho(x),H t (y)} 

To calculate {H (x), %i(y)}, we will first calculate {F H ,G} = {F T ,G} + {F V ,G}. 
Substituting equations (182) and (193) into the bracket {F T , G} yields 

{F T) G} = J d 3 zfV a {g a U T ). (202) 

Assembling equations (186) and (193) into the bracket {F V) G} yields 

{F V ,G} = J d 3 zf°(V a g a )H v + 2 J d 3 zVhf(V c g a )R a c 

+ 2 J d 3 z (V a g a ) V C V7° - 2 J d 3 z Vh (V c g a ) V a V c /° • (203) 

After integrating the last two terms by parts, the identity (V a V c — V c V a )g a = R ac g a 
implies that 

{F v , G} = J d 3 z f {V a g a ) n v + 2j d 3 z VhR a c V c (f°g a ) . (204) 

By integrating the last term by parts, using the identity 2V c R a c = V a R = V a (R — 2 A), 
and recalling that T-L v — V~h(2A — R), the bracket becomes 

{F v , G} = J d 3 z f°V a (g a H v ) ■ (205) 

Combining {F v , G} with {F T , G} and recalling that Ho = Wt + Hv yields 

{F H) G} = J d 3 zf°V a (g a n ). (206) 

Since g a is a three- vector and Ho/Vh is a three-scalar, 

V a (g a n )=d a (g a n ), (207) 
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from which it follows that 

{F H ,G} = J d 3 zf°d a (g a H ). (208) 
To extract the bracket {'Ho(x),'Hi(y)}, first relabel the variable of integration, 

{F H , G} = J d 3 x f°(x)d xa (g a (x)H (x)) , (209) 

then use the identity 

g a (x)H (x) = J d 3 yS 3 (x - y)g a {y)U Q {y) (210) 



to write 



{Ft, G} = J d 3 xd 3 yf(x)g a (y)Ho(y)d x aS 3 (x-y). (211) 
By comparing this expression to (178), we obtain the identity 

{no(x),Hi(y)} = n (y)d x ,5 3 (x - y) . (212) 

{Ho(x),H (y)} 

To calculate {Ho(x),'Ho(y)}, we will first calculate 

{F H , G H } = {F T , G t } + {F T , G v } + {F v , G T } + {F v , G v } . (213) 

It is straightforward to verify that the brackets {F T ,G T } and {F V ,G V } vanish iden- 
tically. To compute {F T ,G V }, substitute equations (182) and (187) into the Poisson 
bracket to obtain 

{F T ,G V } = 2 J d 3 zf-K mn V m V n g° 

- J d 3 z f°g°-L (±HvA + 2VhR mn ir mn ^J . (214) 

Likewise, 

{Fv, Gt} = -2 J d 3 zA-V m V„/° 

+ J d 3 z f°g°^= (\UvA + 2VhR mn ir mn ^J . (215) 
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The sum of the four brackets reduces to 

{F H , G H } = 2j d 3 z (f-K mn V m V n g - g°ir mn V m V n f°) . (216) 

After integrating by parts and recalling that Hi = —2hijV 'k^ k , the bracket becomes 

{F H , G H } = J d 3 z {fWV ig ° - g°WVJ°) . (217) 

Upon expanding the covariant derivatives in terms of partial derivatives and connection 
terms, the connection terms cancel to yield 

{F H , G h } = J d 3 z {fn%g° - g°n%f) . (218) 

To extract the bracket {H (x) : Ho(y)} , relabel integration variables and use the iden- 
tities 

g\x) = J d 3 y5 3 (x-y)g°(y), f{y) = J d'x 5 3 (x - y)f(x) (219) 

to write 

{F H ,G H } = J d 3 xd 3 yf (x)g°(y) (W(x)d x ,5 3 (x - y) - W(y)d yl 5 3 (x - y)) . (220) 
By comparing this expression to (178), we obtain the identity 

{H (x), U (y)} = U\x)d xl 5 3 {x - y) - W(y)d y ,5 3 (x - y) . (221) 

Appendix B: Constraint brackets after imposing % 

We begin with the four constraints After introducing the gauge- fixing constraint 

x = Vh-u(t), (222) 

we need to compute the brackets of each of the five constraints (including x) with X- We 
introduce the smoothing functionals 

F x = J d 3 xf x (x) X (x) , G x ee J d 3 yg x (y) X (y) , (223) 
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where f x and g x are arbitrary rapidly-decaying smoothing functions. We then compute the 
brackets 

{F X ,G X } = J d 3 xd 3 yf x (x)g x (y){ X (x),x(y)} 
{F H ,G X } = J d 3 xd 3 yf(x)g x (y){Ho(x),x(y)} 
{F,G X } = J d*xd*yf(x)gM{ni(x), X (y)}. (224) 



The variation 5F X is 



5F X = J d 3 xf x ^Vhh ij Sh i:j , (225) 



so 



SF 1 i- SF 

°** -f ^ h mn ^£*_=0. (226) 



Likewise, 



Sh mn JX 2 ' 5n 



^2L = g x \Vhh mn , = . (227) 



It follows at once that 

{F x ,G x } = 0. (228) 
Comparing with (224), we obtain the identity 

{X(x),x(y)} = 0. (229) 
We now turn to the brackets of x with the H^. 

• {'Ho{x),x{y)} 

We split {F H , G x } into {F H , G x } = {F T , G x } + {F v , G x }. Assembling equations (182) 



and (226) into the Poisson bracket {F T , G x } yields 



{Ft,G x } = J d 3 zf°g x ^ k k 



(230) 



The bracket {-FV, G x } vanishes identically, so 

{F H ,G X } = J d 3 xf°g x ^ k k . (231) 
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To extract the bracket {H (x), x(y)}, use the identity 

g x (x) = J d 3 yg x (y)5 3 (x-y), (232) 

which yields 

{F H ,G X } = J d 3 xd 3 yf{x)g x {y) l -ir\{x)5 3 {x-y). (233) 
Comparing to (224), we obtain the identity 

{-Ho(x), X (y)} = ^ k k (x)6'(x-y). (234) 

• {n t (x), x (y)} 

From equation (192), it follows that 

{F, G x } = - J d 3 z g x VhVJ l . (235) 
Integrating by parts and using the fact that g x is a scalar, this bracket becomes 

{F,G X } = J d 3 xfVhd i9x . (236) 
To extract the bracket {Hi(x) : x(y)}, use the identity 

g x (x) = J d 3 yg x (y)S 3 (x-y) (237) 

to write 

{F,G X } = J d 3 xd 3 yf\x)g x (y)^h(x)d xi S 3 (x-y). (238) 
Comparing to (224), we obtain the identity 

{Hi{x), X (y)} = Vh(x)d xi 5 3 (x - y) . (239) 

Appendix C: Conformal Decomposition 

Consider a metric g^ u in a number of dimensions d. Denote the determinant of g^ v by g. 
Define the positive conformal factor 

ft=|#| 1/c! >0 (240) 



48 



and the metric 



~g„u = i<?r v v ( 241 ) 

so that 

= Vg^ u . (242) 

By construction, the signature of g^ v is the same as that of g^. Denote the determinant of 
9fiv by g. From the definition of g^, it follows that g = g/\g\, so g = ±1, depending on the 
signature of g^ v . We therefore call g^ u a unit-determinant metric. 

The inverse metrics are related by g^ v = g^Q' 1 . We denote the covariant derivative 
with respect to g^ u by V M , and the covariant derivative with respect to g^ u by V M . The 
connection Y^ v defined by g^ v is 

rj, = \g Xa (ptf va + ^ - d a g^) , (243) 



while the connection defined by g^ is 



1 



f %, = (d»~9™ + du~g,a - d a ~ 9llv ) . (244) 
The connection T* v obeys F* v = F* v — C^ u , where 

= (8% - \g Xa ~9,)j do log n . (245) 
For convenience, we can write Q in terms of a scalar field ip and a constant Q as 

n = Q e 2ip , (246) 

in which case 

= 8fi v <p + #V„¥> - g^V V • (247) 

The Riemann tensor of g^ u is 

R k^v = ^ti^Kv ~~ dv^Kn + r mo t \ v — r^r^ , (248) 

while the Riemann tensor of g^ u is 

R kiiv = d^F KU — d v F Kfl + F ^ KV — F v Jt Rfl . (249) 
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Using = f x + C x v , the Riemann tensor R x Kfll/ can be rewritten as 

pA pA , ria /~t\ ficr _|_ o _|_ -pA /^rcr -per s~t\ 



9uC^ - T x a C^ + Tl K C x a . (250) 



Using 



V/xC^, - V V C X — d^C x v + T x C° KV - T° K C* 



KfJ, ^/J-^KV ' jj,a^ KV flK (TV 

i\ pA /-icr , per /^A 



~ di>C x - T x a C° + T a VK C x , (251) 



it! becomes 



-^ A «^ — ^k/uv + C^Ckv ~ ^ua^K/x + ^t^Kv ~~ ^fC^ ■ (252) 

The Ricci tensor of is i?^ = R x fl \ u ', the Ricci tensor of ^ is R^ v = R x ^\ u - Tracing 
equation (252) appropriately yields 

R^u = Rfiu + C x a C° u - C x a C° x + V\C X U - V V C X X . (253) 

We now express R^ in terms of R^ u and derivatives of <p. Recalling that 5^ = d, we find 

Cl = dV a <p 

= (d + 2)(V^)(V^) - 2^(V^)(V» , (254) 

so 

R»u = R»u + (d- 2)(V^)(V^) - (d - 2)^(V^)(V» 

- (d - 2) VfivV - ^V CT V<V • (255) 

The Ricci scalar for g^ u is R = g^R^; the Ricci scalar for g^ u is R = g^R^. In terms of 
covariant derivatives of ip, we have 

SIR = R - (d - \){d - 2)(V a if)(V a if) - 2(d - l)V a V a (p . (256) 

In three dimensions, the Weyl tensor vanishes, so the Riemann tensor is completely 
determined by the Ricci tensor and the metric via 

Rlkrnn ^ ^ \9lmRkn QlnRkm 9kmR\n QknRliri) 

(9lm9kn- 9ln9km)R- (257) 



(d- l)(d-2) 
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In this case, it suffices to compute the Ricci tensor. When d — 3, our previous formulas 
reduce to 



QR = R — 2( V fc y?) ( V V) - 4V fe V fe v? • (258) 
The condition u ~ uj(t) amounts to ip ~ <£>(i), so i n cosmological gauge we have 

~ ^ , ttR~ R. (259) 
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